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Abstract 

The alternative approach to QCD analysis of the photon structure function F^ is 
presented. It differs from the conventional one by the presence of the terms which in 
conventional approach appear in higher orders. We show that this difference concerns 
also the photonic parton distribution functions. In the alternative approach, the com- 
plete LO analysis of F^ can be performed as all required quantities are known. At 
the NLO, however, one of the coefficient function is so far not available and thus only 
the photonic parton distribution function can be computed and compared to those of 
standard approach. We discuss the numerical difference of these approaches at the 
LO and the NLO approximation and show that in case of i 7 ^ this difference is non- 
negligible and may play an important role in the analysis on photon data of the future 
experiments. 
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1 Introduction 



The concept of the structure of the photon has appeared more then 30 years ago and is 
still very important phenomenological tool for analyzing processes involving photons. It is 
based on the fact that photon fluctuate into various states consisting of leptons, quarks, W ± 
bosons etc. Through the primary splitting 

7 -> qq (i) 

followed by subsequent interactions of the qq pair, the hadronic structure of the photon in 
the perturbative QCD can be revealed. 

This structure can be described, similarly as for the case of nucleon, in terms of parton 
distribution functions (PDF), satisfying DGLAP evolution equations. However, for the pho- 
ton these equations differ from those of nucleon one by the presence of inhomogeneous terms 
which are the direct consequence of ([I]). This inhomogeneity has important implications for 
the QCD analysis of F^. In particular, one has to be careful in the interpretation of the 
order of QCD approximation of physical quantities, like F%, involving PDF of the photon. 
The way of handling this problem in most of analyses of structure function F% , which we call 
the conventional approach, has been critically discussed in p] and the alternative one was 
proposed^. According to pQ, the main source of the confusion concerning the QCD analysis 
of F2 stems from the definition of what the leading order (LO) and next-to-leading order 
(NLO) means in the context of photonic interactions. As a result, the conventional approach 
mixes purely QED effects with those of QCD origin. The alternative approach, compared 
to the conventional one, includes at each fixed order of QCD a few additional quantities. In 
LO the complete analysis in this approach can be performed because all required quantities 
are known. In NLO all but one of these quantities are known but even the remaining can in 
principle be derived from the known calculations [3J. The work on this is in progress [I]. 

In the present article we extend the alternative approach to the case of full F% and 
present detailed numerical results. In particular, we emphasize the numerical importance 
of all additional terms as well as the overall difference between these two approaches and 
show that this difference is phenomenologically relevant and comparable with errors of the 
existing data. Finally we present the global analysis of F^ based on the model of ref. [5]. 

The paper is organized as follows. In the next Section the basic facts concerning PDF 
of the photon are recalled, followed in Section 3 by discussion of the QCD analysis of the 
structure function F^ at the LO and NLO. In Section 4 numerical comparison of the con- 
ventional and the alternative approaches to the QCD analysis of F^ at the LO is performed, 
followed in Section 5 by results of the global fits in FFN Sqjkl model and the conclusion 
in Section 6. 



1 The alternative approach was also adopted for instance in [2]. 
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2 Photonic parton distributions and their evolution 



Similarly as for nucleous, we can introduce the so-called "dressed" photonic PDFs, which 
result from the resummation of multiparton collinear emission from the corresponding "bare" 
parton distribution. As a result, photonic PDFs acquire dependence on the factorization 
scale M. This dependence is determined by the system of coupled inhomogeneous evolution 
equations 



dG^(x,M) 



k G + P Gq + P GG ® G\ (2) 



dlnM 2 

dgl s (x,M) 

dlnM 2 - °NSKq + Pqq ® Q N S 



where 



5 NS = 6n f (<e 4 > - <e 2 ) 2 ) , 5 E = 6n, <e 2 ) (3) 

and the inhomogeneous terms k q a k G result from ([T]) and QCD corrections to it. 
The non-singlet and the singlet distribution functions are defined as: 

W(x, M) = ;>>J(x, M) + q}(x, M)l 
f 

ql s (x,M) = J> 2 - (e 2 ))(q](x,M) + qj(x,M)) 
f 

where / runs over all active massless quark flavors. 

To order a the photon-parton splitting functions k(x, M) and the purely hadronic split- 
ting functions P(x, M) in Eq. (J3J) are given as power expansions in a s (M) 



a 



k q (x,M) = — 



mx) + *m k u {x) + ( *M) 2 m x) + . . . 



2tt 



a 



k G (x,M) = — 
Ztt 



a s {M) 
2vr 



k ( G \x) + 



a s (M) 



2tt 



2tt 



k%\x) + 



Pij{x,M) 



where the splitting functions k g for i — 0, 1, 2 are given in Mellin moments in Appendix A. 
The coupling constant a s satisfies the RG equation: 



(4) 
(5) 
(6) 



— = ftaM) = -- 



(7) 



with fa = 11 - 2//3 and ft = 102 - 38//3. To the NLO, the solution of the equation 
can be approximated in the form: 

Q 8 (aQ ^ 1 _ gi lnln/x 2 /A 2 
4vr ~/3 ln/x 2 /A 2 /3 3 (ln/i 2 /A 2 ) 2 ' U 

Equation ((7j) has infinite number of solutions which can be parameterized by the parameter 
A. The choice of A then specifies a particular solution of ((7|). 

The general solution of the evolution equations (El) can be written as the sum 

P(x, M) = P' had (x, M) + P' PL (x, M) (9) 

of general solution of the corresponding homogeneous equations, called hadronic, and par- 
ticular solution of the full inhomogeneous equations for which we can take the so-called 
pointlike solution which results from the resumation of contributions to the pure QED cou- 
pling 7 — > qq where one gluon, two gluons etc. are emitted. In the non-singlet case the 
diagrams like those in Fig. [1] are ressummed: 



qYs b (x,M a ,M) = 

^ j mi r Jx v \yj J Mi r i ^ 2tt j M 2 t 2 

f 1( fyp(o)(A f 1 ^. P (o) fy\ f M2 dT i a s{ri) 

L V qq \y) Jy W « \w) J ME n 27T 

x r^^M^(o) H r 

J Ml T 2 2tt 2tt J M 2 r 3 




Figure 1: Diagrams defining the pointlike parts of non-singlet quark distribution function. 

The evolution equations (131) can not be solved analytically in x-space and one has to use 
numerical methods. One of them is based on the Mellin transform: 

f(n,M)= [ dxx^ffaM), (11) 
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which converts the equations ([2D into much simpler ones: 
d£T(n,M) 



dlnM 2 

dG 7 (ra,M) 

din M 2 

d<& s (n,M) 
din M 2 



5zk q (n, M) + P qq (n, M)£>, M) + P 9G (n, M)G 7 (n, M), (12) 
fc G (n, M) + P Gq (n, M)E>, M) + P GG (n, M)G J (n, M), (13) 
5 NS k q (n, M) + P qq (n, M)g^ 5 (n, M). (14) 



Evolution equations (|12[) - (|14p are solvable analytically. In the non-singlet case, taking into 
account only k^ and P^ we get [6] : 



Q%s L (n,M) 



Air 



q$*(n, M) 



a s (M) 



a s (M )J 



a 



l-2P^ ) (n)//3 



a 



a s (M ) ) 

qi> had (n,M ) 



a, a- 



n 



-2P^(n)/f3 



Adding also k^ 1 ' and P^ we find: 

47T 



(n,M) 



a s (M) 



J- H - ^iVS 



2^ 



1 - 



a 



:(M) \ 



l-2P^ ) (n)//3 



a s {M )J 



+ 



1 - 



a 



(M) \ 



-2P 9 ( 9 0) (n)//3 



a s {M )J 



a 



a s (M) a s (M ) 



2tt 



2tt 
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NS 



k^(n)-^°\n)-U NS k^(n 
(M) \ 



O.s 



2P^\n)/(3 



(15) 



(16) 



a s (M ); 



~<> had (n, M 



where a s is now given by equation (JTj), where we keep terms proportional to Po and Pi and 



In the flavor singlet case the above formulae become more complicated because 

p(0) P (i) 

and consequently also are the 2x2 matrices. For more details see Appendix B. 

Next step is to transform the distribution function from the Mellin moments back to 
the x-space. This inverse Mellin transformation can be done by integration in the complex 
plane: 



£+ioo 

F(n)x~ n dn 

£— ioo 
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where the integration path lies to the right of the right-most singularity of F(n). 
If the function f(x) is real the above formula may be simplified as: 



7T 



f( x ) = - F(£ + ir])x^~ ir > dr). (19) 



o 



To solve DGLAP evolution equations requires the specification of the boundary conditions 
which contain information on x-dependence of the photonic PDFs at some arbitrary initial 
scale Qq. This boundary conditions are usually parameterized by a set of parameters which 
are determined by the comparison with experimental data. Another parameter which have 
to be chosen is a real parameter £. As is seen from ( Fl8l) . it determines the position of the 
integral path on real axis. In fact, such a choice is arbitrary if £ lies right to the most- right 
singularity and the integral path runs upwards to infinity, but this is no more true if we solve 
this integral numerically. In such a case the optimal shape and position together with the 
upper limit of the integration have to be chosen in order to satisfy the requirements of the 
accurancy or/and computation time. 

3 QCD analysis of the structure function K7 
3.1 Conventional approach at the LO and NLO 

With the explicit solutions for the photonic parton distributions qJj S , S 7 and G 1 at hand 
it is now straightforward to obtain the photon structure function F^x, Q 2 ) given as (in the 
rest of the text we set M = Q) 

^(x,Q 2 ) = q l s (Q)®C q (l) + ^5 NS C,(x,l) + 

<e 2 > £ 7 (Q) <g> C q (l) + <e 2 > G\Q) <g> C G {1) + |- <e 2 > <feC 7 (x, 1) (20) 

where the coefficient functions C q (x), Cg{x) and C 1 entering the F^ in convolution with the 
photon distributions admit the following perturbative expansions 

C q (x,Q/M)= 8(1 -x)+ ^^Cf (*,Q/M) + --- , (21) 

C G (x, Q/M) = ^l C g\x, Q/M) + • • • , (22) 

C,(x,Q/M)= Cf(x,g/M)+ ^4^1c^(x,Q/M) + --- . (23) 

The coefficient functions Cq , and Cj can be found in Appendix A. In order to compute 
the structure function F% to the desired order of QCD one has to provide photonic PDFs 
and coefficient functions to corresponding order of perturbative QCD. 

We find that in the pointlike solution (fT3|) the QCD coupling constant a s appears in the 
denominator. If we take this fact literally i.e. as meaning that the lowest order of pQCD 
pointlike parts of the photonic distribution functions behave like (l/a s ), we conclude that: 
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• at the LO the structure function is related to the distribution functions in the same 
way as in the hadronic case: 

F^ LO (x, Q 2 ) = ql s (x, Q) + (e 2 ) Q) (24) 

where non-singlet and singlet distribution functions satisfy evolution equations includ- 
ing terms k q °^ and only 

• whereas at the NLO 



F ,,NLO { 



x,Q 2 

1 



(e 2 ) 



1+ 



2n~ q 



Cf (1) ® SP>(Q) + ^Cg\l) ® (P(Q) \ + ^Cf 1) 



.(i) 



2tt 



2tt 



where non-singlet, singlet and gluon distribution functions, satisfy now the evolution 
equations including also the splitting functions kq and k^. 

Let us, however, ask the question what happens to the pointlike solution (jToT) and (fTBl (and 
consequently to F^) if we switch off QCD by sending A — > 0. For this purpose, it is sufficient 
consider the non-singlet pointlike distribution function. It is obvious that provided Mq is 
kept fixed when A — > the sum (TTU1) approaches its first term: 



-y,PL 

Ins 



a l(0) 
271^ 



, M 2 
ln Mf 



(25) 



corresponding to the first diagram in Fig. [T] and describing purely QED splitting 7 — > qq. 
This is obvious also from the fact that switching off QCD, keeping Mo fixed, we discard in 
Fig. [T]all diagrams with QCD verticies and are then left with the first one, purely QED one, 
giving (j25p . This shows that the claim q^g oc l/a s is misleading. In fact, this claim can be 
used only as a shorthand for the specification of large Q 2 behaviour of q]jg L ■ Morover, we 
see that in the conventional LO approximation to i 7 ^ (1241) . the effects of QED are partially 



ignored as the splitting function k g ^ is included, whereas C\ u> is not, despite the fact that 
both come from the first diagram in Fig. [U 



(o) 



3.2 Alternative approach 

The alternative approach p] is based on the procedure of separating the QED contributions 
and defining terms "leading" and "next-to- leading" etc. uniquely as terms of the QCD 
origin, starting at a s . 

The crucial role in this approach plays the consideration that q PL is of the order 0(a) 
rather then 0(1/ a s ). Consequently g 7,PL can be written as the sum of two terms: 

q^ PL (x, M, M ) = q^(x, M, M ) + q^(x, M, M ) (26) 
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where the purely QED contribution proportional to a is given as (for single quark flavour 
with unit electric charge) 

qlH(x, M, M ) = |-3 (* 2 + (1 - xf) In ^ (27) 

where M is a free parameter, on which the separation of hadron and pointlike parts in (JH]) 
depends, and qg CD , describing genuine QCD effects satisfies the inhomogeneous evolution 
equation_| 



dq&x,M) dq^ PL (x,M) dq^ L D (x } M) _ 
InM 2 InM 2 InM 2 

^S^+dp 1 "®^^-- (28) 

Consequently, the structure function F^ 1 at the LO approximation includes the QED 
terms: 

i^V, Q 2 ) = Q) + ^\x, 1) (29) 

and the QCD contributions proportional to aa s : 

*Zq5£°(*, Q 2 ) = &^ Q) + ® #S(Q) + ~C?i*, !)• (30) 

By adding also hadronic part of the distribution functions and summing over charges of all 
active quarks and antiquarks, this to full F2 at the LO approximation: 

x r 2,alter\ dj i V ) ~ 

<& s (x,Q) + (e 2 > V(x,Q) + ^Cf(l) ® ql s (Q) + (e 2 ) ^Cf(l) ® JT(Q) + 
fe^f(,,l) + <e>^4 0, M)+fe|;gc«M) + 1). (31) 

The expression ( 13 ip thus differs from the expression in conventional approach (j24j) by 

• the appearance of the corresponding contributions of the photonic coefficient functions 

and C'y , 

• the appearance of the convolution of quark coefficient function Cg with qJj S and £ 7 , 

• the fact that k^ s resp. k^ and Kq are included in the evolution equations for qJj S 
resp. £ 7 . 



2 In the following the dependence on Mq will be supressed 

3 Note that it is consistent to use in (pTTjl the full pointlike (p' PL instead of <?q B £ 





Conventional 


Alternative 




LO QCD 
NLO QCD 


kf\ P(°) 

1.(0) P (o) 1.(1) 1.(1) P (i) r (o) r w r m 


, r(°) 

"T rug 5 ^7 '7 ; 

, 1.(2) 1.(2) r (l) 


, ^q 

r (2) r (2) r (2) 



Table 1: The list of all quantities entering to the definition of F^ in alternative and conven- 
tional approach in leading and next-to-leading order. 



3.3 Alternative approach at the NLO 

At the NLO we add in both the splitting and coefficient functions the terms standing by ol 2 s . 

In TabJUwe list all quantities used in the both alternative and conventional approaches at 
the LO and NLO. At present, all quantities except the coefficient function are known. 
This allows us to perform the global analysis of up to the NLO in the conventional 
approach. In the alternative approach, only LO global analysis is possible to perform. 
Nevertheless, we will make some comments on the calculation in the NLO in the alternative 
approach as well. For this purpose, we focus on the non-singlet part of F^. 

In the standard approach the non-singlet pointlike distribution function satisfies in the 
NLO the following equation: 



dql s (x,M) 
dlnM 2 

(2) 

whereas in the alternative approach, we include also the term proportional to k q 

dql s (x,M) _ 
dlnM 2 

s + + (l) 2fc « 2,(l) ) + {^ p « + diT p «) 8 (33) 

We now investigate the numerical effect of such a modification on the solution of fl33|) . In 
Fig. [5]the comparison of three terms of the expansion of the photonic splitting function k q (x) 
proportional to kf\ kq 1 ^ and k q 2 ^ is shown. The relative importance of this contributions 
depends on M 2 . In Fig. Eh and Fig. [2b they are calculated for M 2 = 100 GeV 2 and 
M 2 = 1 GeV 2 . One can expect that adding the term k^ to the evolution equation has only 
small numerical effect in comparison to the corresponding inclusion of the term proportional 
to ki 1} . 

Solving evolutions (1321) and ( 1331) we find that non-singlet pointlike distribution function 
at the NLO in the alternative approach differs from the conventional one by the term: 
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-0.1 



_Lj 



1^= 100 GeV 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



(a) 



(b) 



Figure 2: Comparison of the first three terms of the expansion of the photonic splitting 
function k(x) for M 2 = 100 GeV 2 (a) and M 2 = 1 GeV 2 (b). 



a a (M) akj 2) ak^pA f a^ y 1 -— \ 

1 _2im\\2-K 2 p An 2 p 2 ) y \a s (M )J J 

which indeed represents only a small correction even in the region x close to 1. The com- 
parison of the non-singlet pointlike distribution functions in both approaches at the LO and 
NLO is shown in Fig. [HJ In this plot we can see the effect of adding the term proportional 
to k^ 1 ' in the case of the LO approximation and the term proportional to k^ 2 ' in the case of 
the NLO approximation (this effect is given by the difference of the solid and doted lines of 
the same colour). 

It is obvious that in the conventional approach the difference between the LO and NLO 
approximations of the parton distribution functions is much bigger than in the alternative 
one. In the conventional approach this difference is given by adding and dominantly 
by kg . Let us recall that in the alternative approach k^ enters the definition of the LO 
pointlike distribution. Consequently, the difference between the two approaches is much 
bigger at the LO than in the NLO. More precisely, the pointlike distribution function at the 
LO in the alternative approach is numerically close to the pointlike distribution function 
in the standard conventional approach in the NLO. The small difference between pointlike 
distribution functions in both approaches in the NLO is due to the small numerical effect of 
the term proportional to k^ discussed above. 



-y,PL,NLO _ -y,PL,NLO 
lNS,alter. lNS,conv. 
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Non-singlet pointlike distribution functions LO x NLO 



0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 

8.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 

Figure 3: Comparison of the pointlike distribution functions in conventional and alternative 
approach at the LO and NLO for M 2 = 100 GeV 2 and A = 313 MeV 

The non-singlet part of in the alternative approach has the following form: 

(g) 2 (1) ® qlsiQ) + *ns^ (cf (s, 1) + ^C7«(x, 1) + {^fcf\x, 1)) (34) 

which differs from the conventional definition ( 1251) by the presence of d}\ Cf } and C ( q\ 
entering through the convolution with the photonic distribution functions. 

As already mentioned, coefficient function C\ has not so far been calculated. This alows 
us to perform in the alternative approach only LO global analysis of F^. To point out the 
difference between the conventional and alternative approach at the NLO, it is sufficient to 
restrict our discussion to the pointlike part of the non-singlet photonic distribution functions. 
The graphic representation of all available contributions entering to the formula (l34"j) is shown 
in Fig. |U 

4 Numerical comparison of the conventional and the 
alternative approaches to F% at the LO 

We now compare the two approaches to the QCD analysis of the non-singlet, singlet and full 
pointlike part of F^ at the LO. We concentrate on the numerical consequences of including 
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NS pointlike contributions NLO 



0.003 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



X 

Figure 4: Contributions to NLO QCD expression for the pointlike part of in the alter- 
native approach for Q 2 = 100 GeV 2 



the additional terms in the alternative approach. We discuss numerical effects of all these 
terms separately as well as the overall differences between these two approaches. 

The individual contributions to the non-singlet and the singlet part of -F^' PL in the 
alternative approach are plotted in Fig. and Fig. [6b and compared to the non-singlet 
and the singlet part of -F^' PL in the conventional approach. 

All contributions included in the alternative approach evolve with Q 2 (except the term 
proportional to C^). As seen from Fig. [Ufc and Fig. EH, in both non-singlet and singlet 
cases, the sum of the terms proportional to the coefficient functions Cj and give the 
positive contribution in the region x < 0.7. The region x > 0.7 is dominated by the negative 
contribution of which is strengthened by negative contribution of C^P close to x = 1. 

The contribution of Cq , entering through the convolution with distribution function q^g 
and S 7 give the positive correction close to x — 1. 

Finally, adding the splitting function into the evolution equations represents an 
important positive contribution for their solutions in the region x > 0.65. This is valid for 
both non-singlet and singlet distribution function. 

Putting all contributions together we compare in the Fig. Efe and Fig. [Uf the non-singlet 
and the singlet -F7' PL in the conventional and the alternative approach. In the alternative 
approach, an d ne higher in the region up to x = 0.8. On the contrary, in the 

region close to x — 1, the values of i^jvs an< ^ * n ^ ne alternative approach decrease 

much faster then in the conventional one. This is mainly due to the inclusion of 

In the Fig. [7] we plot the non-singlet and the singlet parts to investigate the behaviour 
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Figure 5: Comparison of the experimental data of F^ with the solutions of the evolution equa- 
tions in the conventional and alternative approach. We used parametrization of FFNScjkl 
model with parameters of the first row of Tab. [2j 

of the full pointlike in both approaches. In Fig. [7H all the additional contributions are 
plotted separately. As seen from Fig. [7b, the sum of terms ~ ~ Cy and ~ Cq ® q 
rapidely decrease in the region close to x = 1. This is compensated by the term ~ k q l \ 
Putting all terms together (Fig. [Tb) allow us to compare the full pointlike F^ in both 
approaches. 

The significance of numerical difference between both approaches becomes clearer if we 
compare these differences with errors of data on In the kinematical regions where the 
errors are reasonably small we conclude that this difference is numerically important as 
illustrated in Fig. Is evident that the numerical difference between both approaches is 
comparable to the accuracy of measurement of F^. 



13 



NS pointlike contributions 



Singlet pointlike contribution^ 
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x<e 2 >5 E (a/2it)C. 
x<e 2 >(a/27i)cJ i 1| ®E 
x<e 2 >5 r (a/2it)(a /2x)C 



.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



NS pointlike contributions | 
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Singlet pointlike contributions | 
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(e) 



(f) 



Figure 6: Contributions to the LO QCD expression for the non-singlet pointlike part of F^ 
(a), (c), (e) and singlet pointlike part of F^ (b), (d), (f) included in alternative approach for 
Q 2 = 100 GeV 2 . 14 



Structure function F^' 1 



Structure function F)/ 1 




Figure 7: Comparison of the full pointlike i 7 ^ at the LO in the conventional and the alter- 
native approach for Q 2 = 100 GeV 2 . 



5 Fits to data of in the alternative approach 

In this section we present an exploratory study on global analysis of the structure function 
F2 at the LO in the alternative approach and compare it to the FFNSc jkl parametrization9 
0. 



4 The FFNSc jkl parametrization was the first of the serie of the analysis of the structure of the photon 
presented by the CJKL group. The othes could be found in [5] 
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5.1 FFNScjkl parametrization 

In [5] the input scale is chosen small, Q 2 = 0.25 GeV 2 and the input hadronic distribution 
are generated via Vector Meson Dominance (VDM) model 

r^(x,QD = J2^f v MD- (35) 

v Jv 

In practice one takes into account the p° meson and the contributions from other mesons 
are accounted via a parameter k which is left as a free parameter 

r Md (x,Ql)=K^r(x,Q 2 ) (36) 

fp 

where f 2 /4:TT ~ 2. The input distribution functions are parameterized as folows: 

x( p (x,Q 2 ) = 0, (37) 
xv p (x,Q 2 ) = N v x a {l-xf, (38) 
xG p (x,Ql) = N g x a (l-xf (39) 

where ( p , v p and G p stands for the sea-quark, valence-quark like and gluon distribution 
function of the p meson and the following constraints are imposed 

i 

2v p (x, Ql)dx = 2, (40) 



o 



x(2v p (x,Ql) + G p (x,Q 2 ))dx = 1. (41) 



o 



The first is related to the number of valence quarks in the p° meson, and the second one repre- 
sents the energy-momentum sum rule. These constraints allow to express the normalization 
factors N v and N g as functions of a and (3: 

" WTWWY (42 > 

aT(2 + a + f3)-T(2 + a + f3){3-T(2 + a + {3) 

g[a,P) " r(/3 + i)r(i + «)« + r(/3 + i)r(i + a)/? + r(/? + i)r(i + «)' 1 ' 

Together with the parameter k this leaves a, j3 and k as the only free parameters to be fixed 
in the fits to the F%. 

The FFNScjkl parametrization assumes three active massless quarks: u, d and s. The 
hadronic part of the non-singlet and the singlet distribution functions can be described by 
means of (I3"9l in the following way: 

qf s ad {x,Q 2 ) = ~v-r> had (x,Q 2 ), (44) 
Z jMd (x,Q 2 ) = 2v^ had (x,Q 2 ) + 6C(x,Q 2 ) (45) 
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Finally, the contribution to F^ due to the massive c and b quarks is added. It is approx- 
imated by the Bethe-Heitler formula valid for W 2 > <Am 2 h : 



^ h (x,Q 2 ) = 3% 4 h w(x,Q 2 ) 



(46) 



with 



w(x,Q 2 ) = (3 



— 1 + 8x(l — x) — x(l — x) 



Ami 



+ ln( 



1 + 13 , 



A-TJlP 1 

x 2 + (l-x) 2 + x(l-3x) — £ 



x 



Q 4 







[l-x)Q 2 



1 - 



4m| 
W 2 



(47) 
(48) 



5.2 Global fit at the LO 

In this section we briefly present the results of the QCD analysis of the structure function F^ 
in the alternative approach at the LO. We adopt the Fixed-Flavour-Number Scheme model 
of [5] with the same set of the experimental data. 

We have performed minimum \ 2 fit using MINUIT [7J. The set of 182 experimental data 
was used [8]. Systematic and statistical errors were added in quadrature. 

For rif = 3 we took the QCD scale A® = 314 MeV which corresponds to the value 
A^ 4 ) = 280 [S] on the assumption that a s (nY nf ^ = a s (nY rif ~ 1 \ Masses of the heavy quarks 
are taken to be m h = 1.3 GeV and m c = 4.3 GeV 



Approach 


x 2 


K 


a 


P 


Coventional 


357 


1.726 


0.465 


0.127 


Alternative 


591 


0.876 


1.878 


4.788 



Table 2: The \ 2 f° r 182 points and parameters of the fit in conventional and alternative 
approach in the FFNSqjkl model. 

Three parameters a, (3 and k describing the input hadronic distribution function and 
value of \ 2 in both approaches are presented in TabJ2J In the alternative approach in the 
FFNSqjkl we have obtained higher value of \ 2 then in the conventional one [S] and notably 
different set of parameters of the input distribution function. In order to improve the quality 
of the fit one would have to look for the new parametrization. 

Graphical ilustration of the result of the fit is in Fig. [HI For full set of figures, consult 

ma. 
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6 Summary and conclusions 



The basic idea of the alternative approach to the analysis of the photon structure function 
F2 , based on the separation of pure QED effects from those of genuine QCD origin, is 
recalled. This approach differs from the conventional one by reserving the meaning of the 
concepts of the "leading" and "next-to- leading" order to the genuine QCD effects. 

The LO and NLO approximations for in both the conventional and alternative ap- 
proach are discussed in the detail and essential differences are pointed out. It is argued 
that the inhomogeneous splitting function k^ should be included in the evolution equations 
for quark distribution functions of the photon already at the LO. Similarly, the photonic 
coefficient functions Cy (which is of pure QED origin) and Cy as well as quark coefficient 
function should be included in the LO formula for F%- 

At the NLO the alternative approach differs from the conventional one by the presence 
of the inhomogeneous splitting functions k q and k G in the evolution equations and by the 

inclusion of the terms proportional to Cy\ Cy 2 \ Cq 2 ^ and Cq in the formula for F^. Except 

(2) (2) 
for Cy , all these functions are known, and also Cy is calculable with current methods, 

though not quite straighforwardly Pj and there si thus hope that complete NLO analysis of 

F2 in the alternative approach will soon be possible to perform. 

As several quantities which in the conventional approach enter i 7 ^ only at the NLO, 
appear in the alternative approach already at the LO, we compare the latter with the point- 
like parts of the LO as well as NLO approximations in the conventional approach. We 
concentrate on the pointlike parts as this is where the two approaches differ. 

As far as the quark distribution functions are concerned, the LO approximation in the 
alternative approach is very close to the NLO approximation of the conventional approach. 
This reflects the fact that the inclusion of the effect of including the inhomogeneous splitting 
function k^ is bigger than that of the homogeneous splitting functions On the other 

hand the LO approximation in the conventional approach is significantly lower than LO one 
of the alternative approach, in particular for large x. At the NLO, however, this difference 
becomes negligible. The LO approximation in the alternative approach thus captures all 
essential contributions to F^ and in view of the accuracy of the existing data should thus be 
sufficient for its QCD analysis. 
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Figure 8: Fits of experimental data of calculated in alternative approach in the 
FFNScjkl model. 
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APPENDIX A 

Coefficient and splitting functions 

Let us present the formulas of the splitting and coefficient functions used in our calculations. 
They follow the expansions ©-©, ([21])- (El- 

The quark coefficient functions were calculated in MS factorization scheme. We present 
them in x-space. 

The quark coefficient function reads [TT] : 



ln(l — x) 



1 — x 



-3 



1 — x 



4 hire 
1 — x 



+ 6 + 4x 



9 + — )S 1 - x - 2 1 + x) In 

3 x 

(A-l) 



The quark coefficient function Cq 2 '' can be written approximately as [T2]: 
1 / 128 184 

Cf\x, 1) = - — L> 3 - -^-£>2 - 31.1052.Di + 188.641D - 338.5135(1 - x) - 17.74L?+ 
4 V 9 3 

72.24L 2 - 628.8Li - 181 - 806.7a; + 0.719xL^ + L Li(37.75L - 147.1Li) - 

80 / 16 232 

28.384L - 20.70L 2 , - — L 3 + n f —D 2 - —D x + 6.34888A) + 46.85315(1 - x)- 

27 y 9 27 



16 



20 



1.5L| + 24.87L! - 7.8109 - 17.82a; - 12.97x z - 0.185xLg + 8.113L Li + — L + — L z 



9 



where 



Lq = In a:, Li = ln(l — x), D k 



ln fc (l - x) 
1 — x 



(A-2) 



The photonic coefficient function reads: 

C(°\x,l) = (1 - 2x + 2x 2 ) In 



1 — x 



X 



1 + 8ar(l - x) 



(A-3) 



The photonic coefficient function Cy can be represented by the approximate formula: 

C«( x ) = — [(26.67 - 317.5(1 - x))L\ - 72.00L 2 - (1287s" 1 - 908.6)Li 

+ 1598LiL -13.86Lg - 27.74L 2 - 67.33L - 1576 - 2727a; + 3823a; 2 ] (A-4) 

The photonic coefficient functions are given in Mellin moments. The leading and the next- 
to-leading terms are given by [llj: 
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kP s (n) = 3n f ((e 4 ) - (e 2 ) 2 ) Kl , fcgj, (n) = 3n f «e 4 > - (e 2 ) 2 ) k 2 



(i) 



(A-5) 



«i = 2 



n + n + 2 
n(n + l)(n + 2) 



(A-6) 



4 A^2/ \ n t \ 5 no n 2 + n + 2 
« 2 = 3 ((*(») " &(») + 2)2 n(n+1)(n + 2) - 

4 _ . , lln 4 + 26n 3 + 15n 2 + 8n + 4 
— 6i(n) + - 



n 



n 3 (n + l) 3 (n + 2) 



(A-7) 



The non-singlet NNLO photonic coefficient function k NS (n) was obtained by the Mellin 
transform of the h^six) presented in [13] : 



3n / ((e 4 )-(e 2 ) 2 ) 



128 Sf (n) + 6S 2 (n) Sf (n) + 8 S 3 (n) Si (n) + 3S 2 2 (n) + 6 S 4 



n 



n 



112 S f (n) + 3 ft (n) S 2 (n) + 2 S 3 (n) + ^ g Sf (n) + S 2 (n) _ ^ S_i (n) 



9 n 
1353 1262 449.2 

+ 7 7T + 



1445 325.4 



n 



(n + 1) (n + 2) (n + 3) 



n 



Cs + -- 



n n 
C 2 Si (n) S 2 (n) 



- S 3 (n 



?7 



rr 



n 



Cs + C2 Si in) - Si [n) S 2 in) - S 3 (n) 

n \ In 

1169 5Q og S 3 (n + l) + 3Si ( n + 1) S 2 (n + 1) + 2 S 3 (n + 1) 



2 n + n 2 + 1 



744.6 403.2 160 512 



+ 



n 3 n 4 + 9n 5 + Uf V 27 



n + 1 

32 S 3 (n) + 3 Si (n) S 2 (n) + 2 S 3 (n) 



n 



S? (n) + S 2 in) „ Si(n) 40.035 114.4 24.86 53.39 
+258.142 + 18.77 -^-^ + 



17.046 , 

+ Cs + - 

n \ n 



n n 
C2 Si (n) S 2 (n) 



n+1 n+2 n+3 



S 3 (n) + 



rr 



n 



538.8 7 S 2 (n) + S 2 (n) A 

Ca + C2 Si (n) — Si (n) S 2 (n) - S 3 (n) J + 

Si 2 (n + 1) + S 2 (n + 1) 21.55 21.984 64 



n 



26.63 



2 n + n 2 + 1 



270.0 



n + 1 



n 



> + 9n 4 



(A-8) 



where 
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Si(n) = 7£ + ^(n + l), 7^ = 0.577216, 

S 2 (n) = ( 2 -^( n + l), ( 2 = ^ 

o 

S 3 {n) = C 3 + ^"(n + l), Cs = 1.202057 
where ip(n), ip'{n) and ip"{n) are the ^-function and its first and second derivatives. 
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APPENDIX B 

Solution of evolution equations in the singlet case 

To find solution in the case of the evolution equations for singlet distribution function 
S(x, Q 2 ) and gluon distribution function G(x, Q 2 ) is less trivial as these equations are cou- 
pled. It is useful to introduce 2x2 matrices PW anc l 2-component vectors k^: 



p(i) = 



p(*) p(0 
p(*) pO) 



jfeW 



l.0 

k {i) 
G 



(B-l) 



where P ? ( ?, pg\ P s ( ?, PjJ and kf and fc^ are standard parton-parton and photon splitting 
functions. This allow us to write the matrix evolution equation for 2-component vector: 



ql(x,M) 



lP(x,M) 



(B-2) 



in the following form: 



k(x,M) + P(M) ®ql(M). 



(B-3) 



dlnM 2 

If we keep the terms up to P^ and k^ 1 ', the pointlike solution of ( 1B-3I) can be written as: 



q"s PL (n, M) 



47T 



a s (M) 



aJM) 
1 + „ ' U 



2n 



a 



l-2P(°)(n)//3 



a s (M ); 



+ 



-2P<°>(n)//3b' 



a s (M )J 



1 a 
-P(°)(n)27r 



fcW(n)- Ajfe (0) (n)-#Jfe (0) (n) 
2p 



(B-4) 



where 



1 a 
l-^)(n)2^ 



n. 



and U stands for the 2x2 matrix: 



(B-5) 



U 



AfP RP . /> ///> ) . '' m ' , '' 

PO A + — A_ — gPo A_ — A + — 



(B-6) 



where 



p = p(d-Ap (0) , 

2po 



(B-7) 
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p(0) , p(0) , ( p (0) _ p (o)\ 2 4p (o) p (0) 

-T qq T igg =■= y I -fg<? -T I T~ ^iqg *gq 

\t = 2 " • < B - 8 > 

= ± { pm - A T i) 

± A n_ A n V ; 

Hadronic solution is given by: 

2 p(°)(n) 



«r>, M ) = ^ 1 + ( 2^r ~ 2k ) u ) {iM ) qr (n ' Mo) - (B " 10) 
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